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Abstract 

We make a generalization of the type C monomial space of a single variable, which 
was introduced in the construction of type C M-fold supersymmetry, to several vari- 
ables. Then, we construct the most general quasi-solvable second-order operators 
preserving this multivariate type C space. These operators of several variables are 
characterized by the fact that two different polynomial type solutions are available. 
In particular, we investigate and classify all the possible Schrodinger operators re- 
alized as a subclass of this family. It turns out that the rational, hyperbolic, and 
trigonometric Calogero-Sutherland models as well as some particular type of the 
elliptic Inozemtsev system, all associated with the BCm root system, fall within 
the class. 
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1 Introduction 



Finding exact solutions for an equation describing a physical system is one of 
the most important issues in mathematical physics. For quantum mechanical 
systems, the discovery of the underlying Lie-algebraic structure in Ref. [1] has 
stimulated systematic construction of (quasi-) solvable models. The first stage 
of the generalization of the idea to quantum many-body systems was however 
confronted with the difficulty due to the fact that second-order differential 
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operators are in general not equivalent to Schrodinger operators, see, e.g., 
Refs. [2,3,4,5]. The next stage began with the success [6] of an algebraization 
of the known exactly solvable Calogero-Sutherland (CS) models [7,8]. The key 
ingredient in Ref. [6] is to construct the generators of 5l(M + 1) Lie algebra in 
terms of the elementary symmetric polynomials which reflect the permutation 
symmetry of the CS models. The similar methods were employed to show the 
(quasi-) solvability of the variety of the CS models [9,10,11,12]. Finally, the 
complete classification of the models constructed from the 5l{M + 1) alge- 
braization was given in Ref. [13]. It turns out that all the obtained models are 
of the Inozemtsev type, which was shown to be classically integrable [14,15,16]. 
One of the key observations in Ref. [13] is the significance of the underlying 
symmetry of the solvable sector, namely, the subspace preserved by the Hamil- 
tonian. Another aspect of the significant role of the symmetry in the one-body 
case was reported in Ref. [17] in connection with A/'-fold supersymmetry. 

In a previous paper [18], we developed a systematic algorithm to construct an 
A/'-fold supersymmetric system starting from a function space preserved by 
one of the pair Hamiltonians. We then applied it to monomial type spaces to 
find out a new family of A/'-fold supersymmetry called type C. As a byproduct, 
we found that for all finite values of the spectral parameter the characteristic 
equation of the Hamiltonians with type C A/'-fold supersymmetry possesses 
two linearly independent series solutions around the origin in a suitable vari- 
able (cf. Section 5 in Ref. [18]). In this respect, there is a famous mathematical 
theorem, called Fuchs's theorem, which tells a sufficient condition for second- 
order linear differential operators of a single variable to acquire two linearly 
independent series solutions. The applicability of the theorem however is rela- 
tively limited and a generalization to partial differential equations is in general 
difficult. Therefore, the procedure presented in Ref. [18] has a potential pos- 
sibility to be quite efficient in investigating what kind of linear differential 
operators, especially in several variables, can have two or more linearly inde- 
pendent series or polynomial solutions. In this context, it is natural to inquire 
in particular whether a similar approach could be applied to investigate what 
types of quantum many-body Hamiltonians admit the above property. 

In this article, we will show that we can indeed systematically construct a 
family of quantum many-body systems with this characteristic by a suitable 
generalization of the one-body case in Ref. [18]. A key ingredient is how to 
generalize the so called type C monomial space of a single variable, which is 
preserved by the type C A/'-fold supersymmetric Hamiltonians, to several vari- 
ables. This is actually non-trivial problem because there are so many possible 
generalizations and we can hardly know a priori which one is preferable to 
others. For this purpose, we employ significant properties of the type C space 
of a single variable discussed in Ref. [18] as criteria. It turns out that the 
criteria really work advantageously in several contexts. In particular, we can 
mostly avoid the problem caused from the fact that second-order differential 
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operators are in general not equivalent to Schrodinger operators, and we can 
completely keep the same classification scheme as that in the single-variable 
case. 

The article is organized as follows. In the next section, we summarize some 
definitions related to quasi- solvability in order to avoid ambiguity. In Sec- 
tion 3, we show how to generalize the type C monomial space of a single 
variable to several variables. Then in Section 4 we construct the most general 
second-order linear differential operator which preserves the generalized type 
C space. The obtained operators are however not equivalent to a Schrodinger 
operator in general. Thus we present in Section 5 a systematic procedure to 
extract the operators which are equivalent to a Schrodinger operator, namely, 
gauged Hamiltonians from the previously obtained family of the operators. In 
Section 6, we clarify how the gauged Hamiltonians obtained in Section 5 are 
embedded in the general type C operators by identifying the relation between 
the free parameters in the operators. Utilizing the form-invariance under the 
projective transformations, we fully classify the possible Hamiltonians pre- 
serving the generalized type C space in Section 7. The explicit form of the 
potential as well as the solvable sectors are exhibited in each case. In Sec- 
tion 8, we briefly discuss possibility and difficulty in the generalization of the 
type B monomial space to several variables. Finally, we summarize the results 
obtained here and discuss further possible developments of the present work 
in Section 9. Some useful formulas are summarized in Appendix. 



2 Definition 

First of all, we shall give the definition of quasi-solvability and some notions 
of its special cases based on Refs. [4,13]. A linear differential operator H of 
several variables q — {qi, . . . , qu) is said to be quasi-solvable if it preserves a 
finite dimensional functional space V/v whose basis admits an explicit analytic 
form 

HVm^Vm, dimV^ = n(Ar) <oo, Vm = {(f)M, ■ ■ ■ AniN){q))- (2.1) 

An immediate consequence of the above definition of quasi-solvability is that, 
since we can calculate finite dimensional matrix elements Sk,i defined by, 

H<t)k^Y.^kA (A; = l,...,n(AA)), (2.2) 

1=1 

we can diagonalize the operator H and obtain its spectra in the space Va/-, at 
least, algebraically. Furthermore, if the space VV is a subspace of a Hilbert 
space L'^{S) {S C R^) on which the operator H is naturally defined, the 
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calculable spectra and the corresponding vectors in Vj\f give the exact eigen- 
values and eigenfunctions of H, respectively. In this case, the operator H is 
said to be quasi- exactly solvable (on S). Otherwise, the calculable spectra and 
the corresponding vectors in VV only give local solutions of the characteristic 
equation of H. 

A quasi-solvable operator H of several variables is said to be solvable if it 
preserves an infinite flag of finite dimensional functional spaces VV, 

Vi C V2 C • • • C Vaa C • • • , (2.3) 

whose bases admit explicit analytic forms, that is, 

HVm^Vm, dimVA/- = n(A^) < 00, Vu ^ {(t>i{q), ■ ■ ■ ,(t>n{M){q)), (2-4) 

for A/" = 1, 2, 3, . . .. Furthermore, if the sequence of the spaces (2.3) defined on 
5 C satisfies, 

V;:i{S)^L\S) (Af^oo), (2.5) 
the operator H is said to be exactly solvable (on S). 



3 A Generalization of the Type C Monomial Space 

In this article, we shall consider a quantum system of M identical particles on 
a line. The Hamiltonian is thus given by 

1 52 

H = --Y.^ + V{q„...,qM), (3.1) 
^ i=i ^Qi 

where the potential has permutation symmetry: 

V(...,qi,...,qj,...) = V{...,qj,...,qi,...) ^ j ■ (3.2) 

To construct a quasi- solvable operator of the form (3.1), we shall follow the 
three steps after Ref. [13], namely, i) a gauge transformation on the Hamilto- 
nian (3.1): 

H = e^(«)i/e-^(«) , (3.3) 
ii) a change of variables from qi to Zi by a function 2; of a single variable: 

Zi = z{qi) , (3.4) 
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and iii) the introduction of the elementary symmetric polynomials of Zi defined 

by, 



M 

(^k{z) = Zh---^ik {k=l,...,M), (7o = 1 . (3.5) 

ii<---<ik 

Due to the permutation symmetry of the original Hamiltonian (3.1), the 
gauged Hamiltonian (3.3) can be completely expressed in terms of the ele- 
mentary symmetric polynomials (3.5). 

The next task is to choose a vector space to be preserved by the gauged 
Hamiltonian (3.3). The type C monomial space of a single variable z is defined 

by [18], 

V& = (l, ■ ■ ■ , ^\ ^'+^^-') , = M + M , (3.6) 

where A/i and A/2 are positive integers satisfying A/i > A/2 and A is a real 
number with the restriction 

A e M \ {-AA2, -AA2 + 1, . . . , AAi} , (3.7) 

and with A 7^ — A/'2 — l,A/i + 1 if A/i. = 1 or A/'2 = 1. It is however evident 
that there are many possibilities in generalizing the space (3.6) to several 
variables. Hence, we need suitable principles in order to render the possible 
generalizations unique. For this purpose, we shall recall the characteristic fea- 
tures of the type C space (3.6) discussed in Ref. [18]. First of all, the space 
(3.6) decomposes as the direct sum of two spaces: 

V& = vi^^®-^<, (3.8) 
where vj^^ is the type A monomial space of dimension A4 defined by 

Vif,) = (1, ^-^). (3.9) 

Thus, we shall generalize the type C space (3.6) so that a similar decomposition 
is reahzed. The next task is then to find a several-variable counterpart of the 
type A space (3.9). In this respect, we note the fact that the type A space (3.9) 
provides an s[(2) module, which is the foundations of the sl(2) construction 
of quasi-solvable models in Ref. [1]. Therefore, a natural generalization of the 
type A space to several variables is provided by an 5l{M + 1) module given by 

/ ^ \ 

Vif^M = • • : ni e Z>o, < < A4 - 1) . (3.10) 

\ i=l ' 

Indeed, this space is exactly what the quasi-solvable M-body Hamiltonians 
constructed in Ref. [13] from sl{M + l) generators preserve, and reduces to the 
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type A space of a single variable (3.9) when M = 1. From the decomposition 
criterion, the generalized type C space we shall take thus have the following 
form: 

Vi?k;M = Vi;^!M ® fx(z) Vi^U ■ (3.11) 

The function fx in Eq. (3.11) must be fx{z) ~ when M = 1 in order that 
the space (3.11) reduces to Eq. (3.8). Next, we recall the fact that the type C 
space (3.6) is form-invariant under a special projective transformation. More 
explicitly, the space (3.6) satisfies 

In Ref. [13] it was shown that the invariancc under projective transformations 
of the vector space preserved by the Hamiltonians plays crucial roles in con- 
structing quasi-solvable quantum many-body systems. Thus, we shall require 
that the generalized type C space has a similar property to Eq. (3.12). To see 
how the space of the form (3.11) behaves under the special projective trans- 
formation, we first note that the elementary symmetric polynomials (3.5) are 
transformed according to 

ak{z) ^ (Tk{z~^) = aM-k{z) crM{z)~^ . (3.13) 

Prom this formula, we easily see that the type A space of several variables 
defined by Eq. (3.10) has the following invariance property: 

Vi?jMW - = a-J-^^-'^ vifj^W • (3.14) 

Hence, the space (3.11) transforms as follows: 

It is now clear that the space (3.11) has a suitable form-invariance similar to 
Eq. (3.12) if and only if the function fx in Eq. (3.11) satisfies 

^^'-^'h{z-') = fAf,-M.-x(z) . (3.16) 
Actually if this is the case, we have 

V)^k;M[^, A] ^ <7^^-^Vi^^v,;M[^-S A] = V)^k;M[^,-A/'i - J^2 - X] , (3.17) 

which reduces exactly to Eq. (3.12) when M = 1. The solution of Eq. (3.16) 
is easily found if we note that the dependence on J\fi — A2 only comes from 
the variable ctm and take the formula (3.13) into account. Only the possible 
solution for fx, which reduces to z^ when M = 1, is given by 

fx{z)=aM{z)\ (3.18) 



6 



In summary, we have generalized the type C monomial space (3.6) of a single 
variable to several variables by the requirement of the decomposition structure 
(3.11) and of the form-invariance under the special projective transformation 
(3.17), and finally obtain 

"KaA,JV2;M = '^Xi;M ® "Ka^Im ■ (3.19) 



4 Construction of Quasi-solvable Operators 



In this section, we shall construct the general quasi-solvable operators of (at 
most) second-order leaving the type C space (3.19) invariant. For the first- 
order differential operators, we have the following set of independent ones: 



^ = E^ (z = l,...,M-l), (4.1a) 
= Eij (i = 1, . . . , M; J = 1, . . . , M - 1) , (4.1b) 

d 

(7m^ — = Emm ■ (4.1c) 

d(JM 

In the above and hereafter, we employ the convention that the Latin indices 
take integer values from 1 to M while the ones with bar take integer values 
from 1 to M — 1. We note that Eqm and E^m cannot preserve the elements 
a"^ . . . (T^^fT^ inside the space V^^jvi;^^ due to the restriction on A, Eq. (3.7). 
We also note that any first-order operator which raises the degree cannot 
preserve the space V^^j^^.^- 

For the independent second-order differential operators, we have the following 
set: 

E,-,E,j {k > I) , (4.2) 



q2 

^iQ^_ = Efk^oi {k > I) , (4.3a) 

(^MT. 7. — — EmmEqi , (4.3b) 

d f d \ 

c"M-^ ) = Emm,om , (4.3c) 



daM V daM 



- - 

q^Tq^. = ^i'kEji - SfkE,i (i>j;k>l), (4.4a) 
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daMdcTk 



'M 



d 



'da 



M 



d 



da 



M 



Emm Emm — Emm 



-X]=E 



MM,iM 



(4.4b) 
(4.4c) 
(4.4d) 



( ^ d \ ( ^ d \ 

- 1 - g (a + AT, - 1 - g . . (4.5) 



The operators (4.3c), (4.4d), and (4.5) cannot be expressed as polynomials 
of the set of the first-order operators (4.1). In the case of M = 1, only the 
four independent operators, namely, Eqs. (4.1c), (4.3c), (4.4c), and (4.5) with 
i — M — 1 exist, which reproduces the result for the type C quasi-solvable 
operators of a single variable given in Ref. [18]. We also note the important 
fact that all the operators listed above preserve separately both the subspaces 
of '^Mi!j\r2;M Eq. (3.19), as in the case of M = 1 [18]. The most general linear 
differential operator of (at most) second-order which leaves the type C space 
(3.19) invariant is thus given by 

H—— Ao,oo-5'io,oo — X! ^tj.k,uiEij,kEi,i — ^Amm,p,mEmm,p,m 

i H>v,k>l P 

~ X/ ^MM,ij.kE mmE ij^i — Amm,mmEmmEmm 

ix,k 



+ X] ^tikE^k + Emm Emm — Co , 

H,k 



(4.6) 



where the coefficients A, B with indices and Cq are real constants. The Greek 
indices without(with) bar take integer values from to M(M — 1), respectively. 
The summation above is understood to take all the possible integer values 
indicated by the corresponding type of indices. In terms of the variables Ufe, 
the operator H. is expressed as 



M-1 



fcJ=l 



^ = - E [-^0((T)(7fc(7, + Afe;(a) 

+ ^Mfc(0-)J 



92 



M-1 



daida-i 



E [2^o(<7)(7^ 



fe=i 



C^M 



daMdai 



AQ{a)aM + Amm(c^) 



(7m 



M a 

+ E [(-^ + ^ - 3) Ao(<7)<7fe - Bkia)] — - C{a) , (4.7) 



fe=i 



where Aq, A^i, B^, and C are polynomials of several variables of at most 
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second-degree given by 



Ao{a 




(4.8a) 




(4.8b) 



M 




(4.8d) 



(4.8c) 




— (A — I)Amm(c) + (AAmm.mm — Bmm) c"m , 



(4.8e) 



(M - l)(A/'2 + A - l)Ao((7) + Co . 



(4.8f) 



Except for the operators £^io,oo given by Eq. (4.5), all the operators in Eq. (4.6) 
leave invariant both the subspaces in Eq. (3.19) for arbitrary positive-integer 
values of Mi and A2- Therefore, the operator ii is solvable if all the coefficients 
^io,oo vanish, or equivalently, if 



5 Extraction of Schrodinger Operators 

In the preceding section, we have constructed the most general quasi-solvable 
second-order operator if preserving the type C space (3.19). By applying a 
similarity transformation on if, we may obtain a family of quasi-solvable op- 
erators. However, second-order linear differential operators of several variables 
arc in general not equivalent to Schrodinger operators. This fact is one of the 
most obstacles in constructing a quasi-solvable quantum many-body systems. 
Recently in Refs. [13,19], it was shown that the amount of the difficulty can 
be significantly reduced by considering the underlying symmetry of the invari- 
ant space of quasi-solvable operators. In this section, we will show that the 
symmetry consideration is indeed quite efficient in our case too. 

To begin with, let us briefly review the consequence in the case of the type 
A space (3.10) (see Ref. [19] for the details). The symmetry which transforms 
the type A space to itself is GL{2, R) of hnear fractional transformations of Zi 
introduced by 



Ao(a) = 0. 



(4.9) 



azi + (3 



(a, 7, 5 e M ; A = a5 - /?7 7^ 0) . 



(5.1) 
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Then, it is easy to see that the type A space (3.10) is invariant under the 
GL(2,M) transformation induced by (5.1): 



M 



(5.2) 



1=1 



Prom this invariance, we can conchide that any operator which preserves 
the type A space must be shape invariant under the GL{2, M) transformation 
as follows: 



M 



M 



i=l 



1=1 



(5.3) 



where ai stands for all the free parameters involved in Hj^k- particular, 

all the second-order operators Hj^'' which preserve the type A space and are 
equivalent to the Schrodinger operators as well must have the latter shape- 
invariance. From this requirement we can derive that Hj^^ must have the 
following form: 



(5.4) 



where Ck and Rk are constants and ~ ^) shape-invariant 

operators satisfying Eq. (5.3) as the following: 



H 



(0) 

A4 



H 



(1) 



H 



(2) 



M 



52 , A4 - 2 



M 



d 



i=l 



dzi 



(■A4-l)(.A4-2) ^^(,).,^, 

-■-^ i=l 

1 uZn ^ -1 



Eof'(^.), 



i^j Zi Zj dz. 



i=l 



i^j Zi Zj 



i=l 



(5.5) 
(5.6) 



(5.7) 



The functions 0^\z) (n = 2,4) in Eqs. (5.5)-(5.7) are polynomials of (at 
most) nth-degree and are transformed under the GL(2,M) transformation ac- 
cording to 



0^^\z) ^ A-"/2(^^ + S)''Ol^\z) . 



(5.8) 
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Keeping these facts in mind, let us come back to consider the type C case. 
As wc have mentioned previously, all the operators in Eq. (4.6) and hence H 
preserve each of the subspaces in Eq. (3.19) separately. In other words, the 
second-order operator aj}^ ^^^Ha^^ ^^'^ {k = 1, 2) leaves the type A space Vj^^ 
of dimension invariant. On the other hand, the results for the type A case 
just summarized tell us that these operators are equivalent to a Schrodinger 
operator if and only if they have the form given by Eqs. (5.4)-(5.7). Therefore, 
the general type C operator H (4.7) is equivalent to a Schrodinger operator 
of several variables if and only if the following conditions are simultaneously 
satisfied (up to constants): 



(A) 



' M 



{A) 
Ml 



(5.9) 



The exphcit form of the operator ctmH^mI^m is calculated as 



M pa M 

i=l ^^i i=l 



A4-2 + (M-l)c. 



+ 2X0i'\z,)zr'-0','\z^, 



d_ 
dzi 



i^j Zi Zj dzi 



M 

M 

-E 

i=l 



0f\z,) + 2X0r{z^)zr' 



ik). 



Zi Z-i 



12 



-(A4 - 2 + 2(M - l)cfe)0f "(^0 



+ ^'-^^[''-'^'' XOf\z.)zr^ + A(A + l)Of\z.)zr^ 



^0?\z.)-XOf\z.)z-^ 



-Rk. (5.10) 



It is evident that the condition (5.9) is satisfied if and only if the coefficients of 
the differential operators of each order in^^^^ and o-l^Hj^^a]^^ coincide with 
each other (up to a constant for the zeroth-order) . Prom Eq. (5.10) we thus 
obtain the following set of conditions: 



0?\zi)-0','>{z,) 



0^\z, 
(1), 



^0f\zi) = A{zi), 



A'{zi) + 2XA{zi)z-^ , 



M 

E 

i=l L 



0'^\z,) 



(5.11) 
(5.12) 



M 

E 



'M2-M1 

12 



{H -^ + 2{M -l)c)A'{zi) + 



M2-2 + {M -l)c^ 



XA{z,)z.' 
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+A(A + l)A{zi)z- 



M 



+ const., 
(5.13) 



where c = Ci = C2. By substituting Eq. (5.12) for O2 {z.j) in Eq. (5.13), the 
third condition can be expressed solely in terms of 02\zi) and A[zi): 



ETA/^-M (1),, . N -1 
0\ [zi] + \0\ '{zi)Zi 

i=\V ^ 



M 

E 

i=l 



12 



(M- 2^ + 2(71^-1)0) A"(^i) 



+ M-2M + (M-l)c ^^.^^^^^_, ^ ^^^^ _ ^^^^ 



M 



^A'(z,)+2A>1(;.,);.-^ 



2^7' Zq 



+ const. (5.14) 



On the other hand, we note that the l.h.s. of Eq. (5.12) is a polynomial of at 
most second-degree. Hence in the polynomial A{z) of at most fourth-degree 
defined by Eq. (5.11), the coefficients of the fourth- and zeroth-degree must 
vanish in order that Eq. (5.12) is met. We can thus put 

A{z) = a^z^ + a^z^ + axz , 0^2\z) = h^? + h^i^ z + h^o\ (5.15) 

Substituting the ansatz (5.15) in Eq. (5.13), we find the third condition (5.13) 
is fulfilled if and only if 



Ml - 2N2 - 2A - (M - l)c 

O2 = 7, «3 , 



A/l - 2A + (M - l)c 



(5.16a) 
(5.16b) 



Finally from Eqs. (5.9), (5.10), (5.15), and (5.16), we find that the second- 
order operators if which preserve the type C space (3.19) and are equivalent 
to a Schrodinger operator as well must have the following form: 



M 



M 



d 



M 



i=l 



A{zi) d 
Zi Zj dzi 



- C(a(; 



where A{z) and B{z) are given by 

A{z) = a^z^ + a2Z^ + aiz , 



B{z) = -(A/' + A- 3 + 2(M - l)c)a^z^ 

+ (^'i'^ - [M - 2 + (M - l)c\a2)z - (A - l)ai 



(5.17) 

(5.18) 
(5.19) 
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Furthermore, the function C in Eq. (5.17) is calculated as 



M 



- 2 + 2(M - l)c) Y.^"{z,) 



Ml -I 



i=l 



M 



M 



M 



+ Ri 



= (M-l)(AA2 + A-l)a3E^i + co, 



(5.20) 



where the constants Ri and cq are related with each other so that Cq in 

Eq. (5.20) exactly coincides with that in Eq. (4.8f) through the following 
relation: 

Co = ~ ^ (M - 2 - (M - l)c)Ma2 - ^ Mh^^^ + i?i . (5.21) 

In fact, it is easily shown that the gauged Hamiltonian (5.17) can be cast into 
a Schrodinger operator by a gauge transformation 



H = e-^^e^ 

1 M Q2 



M 



2gag,2 + 2g 



dqf 



C(a(; 



(5.22) 



if the gauge potential W is chosen as 

^ r Biz ) 1 ^ I 



M 



= Eln 



Zi Z. 



■3\ ' 



(5.23) 



2A{zi) ' 4 

and the function z{q) which determines the change of variables satisfies 

z'iqf = 2A[z{q)) = 2{a3z{qf + a2z{qf + a,z{q)) . (5.24) 

To make the correspondence between the results here and those in Ref. [18] 
more transparent, let us decompose B{z) as 



(5.25) 



so that the coefficients of the second- and zeroth-degree in Q{z) coincide with 
each other: 



(5.26) 



The parameter f3 in Eq. (5.26) is given by 



(5.27) 



13 



Then, if wc put bi as the coefficient of the first-degree in Q{z) and take M — 1, 
we can easily check that Eqs. (5.17)-(5.26) completely reduce to the corre- 
sponding quantities in Ref. [18] as they should. 

It is apparent from the construction that the Hamiltonian (5.22) preserves the 
space Vj^^VziM defined by 



V 



(5.28) 



Hence, the Hamiltonian H can be (locally) diagonalized in this finite dimen- 
sional space (5.28). We shall thus call the space (5.28) the solvable sector of 
H. To obtain the explicit form of Vj^^^j^^^.^j^, we first arrange the gauge potential 
(5.23) in a more tractable form: 



- M 



— ^\n\zi\ }^ln 



1=1 



+ 



M „ 

/3E/ 

i=i 



4 

M 

c^ln 



A{z,) 



where a new parameter a~ is introduced by 



Zi Z' 



■3\ 1 



(5.29) 



1-2A 



a 



(5.30) 



Then, from Eqs. (3.19), (5.28), and (5.29), we have 



V 



V 



(1) 



Ni;M 



e V 



(2) 



where the new gauge factor e " is defined by 



M 



-U 



i<j 



M 

n 



A{z,) 



i(Ar-l+2(M-l)c) 



X exp 



M „ 

-/3E / 



(5.31) 



(5.32) 



6 Gauged Hamiltonians in the general type C operators 



In the preceding section, we have constructed the type C gauged Hamiltonian 
H (5.17) almost independently of the explicit form of the most general type 
C operator H given by Eq. (4.7). In fact, all that we have employed in the 
previous construction is the known results for the many-body type A gauged 
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Hamiltonians and the fact that H, preserves the two type A spaces in Eq. (3.19) 
separately. In this section, we will thus show how the gauged HamiltonianiT is 
embedded in the general type C operator 7i. For this purpose, we shall express 
H in terms of cxfc. The necessary formulas are summarized in Appendix. 

First from Eqs. (4.8f) and (5.20), we easily see that Aq having the form (4.8a) 
must be 



Ao((7) = 03(71 . 



(6.1) 



With the aid of Eqs. (A.l) and (6.1), the second-order operator in Eq. (5.17) 
is expressed in terms of Uk as 



M q2 M q2 



i=l 



k,l=l 



(6.2) 



where Am are polynomials of second-degree as the following: 



Aki{cr) = - as 



X (/c - m + l)ak-m+icri+m 

k+1 

+ - I - 2m + l)(Tk-m+l'7l+m 

m=2 

k 

k (TkOi + ^{k-l- 2m)(7fe_^(7i+, 

m=l 

k-l 

- oi X (A; - / - 2m - l)ak-m-i(^i+m ■ 

m=0 



+ 02 



(6.3) 



Comparing Eq. (6.2) with the second-order operators in Eq. (4.7), we imme- 
diately obtain 



A-kiW) = AM{<y) + Ark{(7) 

^fefc(o-) = A'fclo") , 
AhdkW)<^M = A'A/lo") +w4ji^(cr) , 

Amm (o") ctm = Amm (o-) • 



{k > I) , 



(6.4) 
(6.5) 
(6.6) 
(6.7) 



The r.h.s. of Eqs. (6.4) and (6.5) indeed have the same forms as Eqs. (4.8b) 
and thus completely determine all the polynomials A^i {k >J). To factor out 
(7m in the r.h.s. of Eq. (6.6), the following formula is useful: 



M-k 



X (M - A; - 2m - l)aM-m-icrk+m ^ {k - M - l)ak-i(TM ■ 



(6.8) 



m=0 



Prom Eqs. (6.3)-(6.8), the first-degree polynomials AMk of the form (4.8c) 
read 



[a) = —2{k + l)a3(7fc_,_i -|- 2ka2(Tk ~ 2(A; — M — l)ai(7fe_i , 



(6.9) 
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AuMicr) = Ma2(JM + oictm-i ■ 



(6.10) 



Similarly, with the aid of the formulas (A. 2), (A. 3) and Eq. (6.1) the first-order 
operators in Eq. (5.17) are expressed in terms of cr^ as 

-2^B{zi)—-2c\^^ 



i=l ijkj 



M Q 

= Y^UAf + X- 3)Ao{a) - B,{a)] — , (6.11) 

k=l 

where Bk are polynomials of first-degree as the following: 

Bk((T) ^(k + 1)(X + X - 3 + kc)a3ak+i 

- ^ (^[AA + 2A - 4 + 2(A; - M)c]a2 + 4/3)afc 

+ {k-M -l)(^X-l + {k- M)c) ai(7jk_i . (6.12) 

In particular, Bm can be arranged with the use of Eq. (6.10) as 

Buia) = -(A - 1) Amm((7) - y ((A^ - 2)a2 + 4/3)(7m ■ (6.13) 

It is now apparent that Bk are all in agreement with the form of Bk given by 
Eqs. (4.8d) and (4.8e) and thus 

Bk{a)=Bk{a). (6.14) 

It should be pointed out that after determining all the parameters A^kji from 
Eqs. (6.4), (6.5), (6.9), and (6.10) we can always adjust the parameters S^j^ and 
Eofc in Eq. (4.8d) and Bmm in Eq. (4.8e) so that Eq. (6.14) holds. Therefore, 
the type C gauged Hamiltonian (5.17) in fact constitutes a subclass of the 
general type C operator (4.7). 

Prom Eq. (6.1) and the discussion following Eq. (4.8f), it is apparent that the 
gauged Hamiltonian H is solvable if 

a3 = 0. (6.15) 

In this case, the polynomial C given by Eq. (5.20) is a constant and hence the 
many-body Hamiltonian (5.22) is of the supersymmetric form up to this irrel- 
evant constant. The gauge potential W can be regarded as a superpotential. 
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7 Classification of the Models 



We shall now exphcitly compute the type C models. The potential term in 
Eq. (5.22) is explicitly calculated in terms of Zi from Eqs. (5.18)-(5.23) as 



M 



1 4:A{z, 



X 



M 



-a3{M,J\fi,X)J2zi 

i=l 



i<j 



[Zi - Zjf 



(7.1) 



where, and in what follows, Vq denotes an arbitrary constant, and the coupling 
constant a^{M,Mi,\) is given by 



a^{M,Ni,X) = - 2M(M + A) + 2(AA+ A)(M - l)c 



as- 



(7.2) 



+ 2M(M - - 2(M - l)c - 1 

From Eq. (5.24), the change of variable is determined by the elliptic integral: 

dz 



H<i-%) = J 



yj2{a3Z^ + a2z'^ + aiz) 



(7.3) 



In contrast to the systems constructed from the si{M + 1) generators, our 
present models do not have full GL{2,M.) invariance; the conditions (5.12) 
and (5.13) apparently break the GL(2,]R) symmetry. However, there remains 
residual symmetry unbroken, namely, the invariance under scale transforma- 
tions, /9 = 7 = in Eq. (5.1), and under special projective transformations, 
a = 6 = and f3 = 7(7^ 0) in Eq. (5.1). The latter invariance is as it should 
be since we have made the generalization so that Eq. (3.17) holds. Due to 
the invariance, we can classify the models completely the same way as in the 
single- variable case presented in Ref. [18]. This is another virtue of our choice 
of the space (3.19) based on the invariance criterion. It is readily shown that 
A{z) can be cast into one of the canonical forms listed in Table 1 by the com- 
bination of the aforementioned scale and projective transformations. For the 
detailed treatment of the eUiptic integral (7.3) in each case, see Ref. [18]. 

Furthermore, we note that from Eq. (5.24), a rescaling of the coefficients Oj, 
/5, Co by an overall constant factor u > has the following effect on the change 
of variable z{q): 



z{q]vai,v(5,vci)) = z{y/u q;ai, P,co) 



(7.4) 
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Tabic 1 

Canonical forms for the polynomial A(z), (5.18). 



Case 


Canonical i^orm 


1 


2z 


2 


±2uz'^ 


3 


±2vz{l -z){l- mz) 


4 


2vz{l — z){m' + mz) 


5 


\vz{z'^ + 2{l - 2m)z + l) 



In this Table, z/>0, 0<m<l and m' = 1 — m (m 7^ 1 in Case 3 and m 7^ 0, 1 in 
Case 4 to avoid duplications). 

From this equation and Eqs. (5.18)-(5.20), (5.22), and (5.23) we easily obtain 
the identities 

W{q ; vai, v(3, vcq) = W{^/i'q ; a^, /3, Cq) (7.5a) 
V{q ; vai, v(5, vcq) = v V{^/v q;ai,(3, cq) . (7.5b) 

We shall therefore set i/ = 1 in the canonical forms in Cases 2-5, the models 
corresponding to an arbitrary value of u following easily from Eqs. (7.4) and 
(7.5). It should also be obvious from Eq. (7.3) that the change of variable 
z{q), and hence the potential V determining each model, are defined up to the 
transformation q ±(g — go)) where go ^ is a constant. In the following, 
another new parameter is introduced by 

a+ ^^f2-^f^ + X + ^. (7.6) 

The condition for solvability (6.15) implies that the models corresponding to 
the first and second canonical forms in Table 1, and to the third one with 
m = 0, are not only quasi-solvable but also solvable. 

As we will see below, the potentials V in all the cases but the irrelevant Case 
2 have the singularities at g^ = and at integer multiples of a real number Q, 
qi = nfl (n G Z), where is a submultiple of the real period of the potential 
y or = 00 when the potential is non-periodic, as well as g^ = ig^ for i ^ j. 
Hence, the Hamiltonians are naturally defined on 

< gM < • • • < gi < ^ • (7.7) 

The potential form of the one-body part in each case is completely identical 
with that of one-body type C A/'-fold supersymmetric potential V~ in Ref. [18]. 
Thus, the normalizability conditions coming from the one-body part, which 
restrict the values of the parameters and in each case, are completely 
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the same as those shown in this reference. Only an additional condition arises 
from the two-body part. The finiteness of the norm of the two-body wave 
function in the solvable sector Vj^^^j^/^-M leads c > —1/2, where c denotes the 
coupling constant of the two-body interaction appeared in the last line of 
Eq. (7.1). 



7.1 Case 1: A{z) = 2z 



Change of variable: z — (f. 



Potential: 



o2 M 
^ 1=1 



a-{a- -1)^ 1 



M 
i=l 



M 



+ c(c-i)E 

i<j 



+ 



+ Vo. 



(7i 



Solvable sectors: 

M 



i<j \ i=l 

M / M 

vZm = n k.^ - ^1 

i<j \ i=l J \ ^ 1=1 



M 



M 



(7.9a) 
(7.9b) 



where V^^lj\^ {k = 1, 2) are defined in Eq. (5.31). This case corresponds to the 
rational BCm type Calogero-Sutherland model [20]. The one-body part of the 
potential is only singular at = and hence a natural choice is Q = oo. 



7.2 Case 2a: A{z) = 2z^ 



Change of variable: z — e^. 



Potential: 



M I 

V = c{c-l)Y. —ry r + Vo . (7.10) 
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Solvable sectors: 



V(2) 



Parameter: 



M / M _ \ 

i<j \ i=l / 

M / M 



n 



n 6=^* -e'* 



.2? 



He 

.1=1 



(l-a--/9)% ■j>(^) 



TV - 1 + 2(M- l)c 



.2? 



(7.11a) 
(7.11b) 



(7.12) 



We note that the potential (7.10) depends on neither the parameter a~ nor (3 
appeared in the solvable sectors (7.11). The dependence on these parameters 
only comes from the middle part of the each solvable sector. These parts 
arc nothing but the wave functions corresponding to the free motion of the 
center-of-mass due to the translational invariance of the potential (7.10) and 
thus irrelevant. The remaining factors in Eqs. (7.11a) and (7.11b) coincide 
with each other and the dependence on is not relevant since this case 
is solvable. Therefore, Case 2 only reproduces the well-known result of the 
(hyperbohc) Am-i Sutherland model [8]. 



7.3 Case 2b: A{z) = -2z'^ 



The formula for the potential for this case can be easily deduced from those of 
the preceding one by applying Eqs. (7.5) with z/ = — 1. Apparently, this case 
only shows the usual solvability of the trigonometric Am-i Sutherland model. 



74 Case 3a: A{z) ^ 2z{l - z){l - mz) 



Change of variable: z — siP q. 



Here, and in the following cases, the Jacobi elliptic functions have modulus 
k = \/m. 



Potential: 



mQ;+(a+ - 1) , «"(«"- 1) 



i=i ^ i=i Qi 
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2 £tcn2g, 



+ c c 



i<j 



+ 



_sn'^{qi-qj) sn'^{qi + qj)_ 



+ Vo. (7.13) 



Parameters: 



:l--UAf-l + 2{M-l)c) 



(7.14) 



Solvable sectors: 



V, 



(1) 



Yl |sn^ qi - sn^ 

i<j 

M 

xl[[{snq,r-{cnq,f-{dnqi)f'' 



sn q 



(7.15a) 



V(2) 



M 



Yl sn^ - sn^ qj 



i<3 



M 



X \\\{snqif~" {cnqiY (dng^)^ 



i=l 

We note that from the formula 

1 



sn^ q 



(7.15b) 



m 
cn^ g 



= — msn g 



+ 



+ 



m 



sn^ g sn^ 2g dn g 



(1 + m), 



(7.16) 



we can regard the above model (7.13) as a deformation of an Olshanetsky- 
Perelomov type potential associated with the BCm root system [20]. In con- 
trast to the ordinary elliptic integrable models, the potential function here is 
given by the Jacobian function instead of the Wcicrstrass p function. We thus 
call the potential (7.13) Jacobi- elliptic BCm Inozemtsev model. The one-body 
part of the potential have real singularities at the points qi — nK{m) {n e Z) 
where K{m) is the complete eUiptic integral of the first kind defined by 



K{m) = / 
Jo 



r^f^ dt 
'0 \/l — m sin^ t 

Hence a natural choice is O = Kim). 



(7.17) 



As previously mentioned, when m = the polynomial A{z) is of second degree 
and hence the Hamiltonian H is solvable. The formulas for the potential and 

the solvable sectors arc obtained from Eqs. (7.13)-(7.15) by setting m = 0, 
m' = 1, and (sng, cng, dug) = (sing, cos g, 1). With the aid of the formula 
(7.16), we immediately see that the potential (7.13) in this case is nothing but 
the trigonometric BCm Calogero-Sutherland model. 
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7.5 Case3h: A{z)^-2z{l-z){l-mz) 



As in Case 2b, the formulas for this case can follow from those of the preceding 
one using Eqs. (7.5) with i/ — —1. The following well-known identities [21] may 
be of help in this case: 



.sn(g;m') 

sn(ig;m)=i^ -, cn(ig;m) = - 

CQ.[q]m') cn[q;m') 



1 A \ dn(5;m') 

^ , dn(ig ; m) - 



cn(g ; m') 



The real singularities of the one-body part of the potential are now at Qi — 
nK{m') {n e Z) and thus we naturally choose as Q = K{m'). The value 
m = yields again solvable model, whose potentials and solvable sectors fol- 
low from Eqs. (7.13)-(7.15) by setting m = 0, m' = 1, and (snig, cnig, dnig) = 
(isinhg, coshg, 1). Clearly, the potential in this case corresponds to the hy- 
perbolic BCm Calogero-Sutherland model. 



7.6 Case 4: A{z) = 2z{l - z){m' + mz] 



Change of variable: z — cn^ q. 



Potential: 



V ^ 



M 



m'a ia — 1) 



i=l 



M 

E 

i=l 



cn^ qi 



+ 



r(/9- 



1).^ 1 



E 



m'/3+(/?^ 



1).^ 1 



i=i sn- qi 

M 



S qi 



+ c(c-l)E 



sn^{qi-qj) sn^{qi + qj) 



+ Vo. (7.18) 



Parameters: 



P^-T/3-l{^f-l + 2{M-l)c). 



(7.19) 



Solvable sectors: 



M 



i=l 



X 



M 

n[(cngO""(sng,)^"(dng,)^ 



cn^ q 



i=l 



"KS;M =]lW(li- sn^ qj 



M 



1=1 



(7.20a) 
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X 



M 



1=1 



cv? q 



(7.20b) 



We again obtain a Jacobi-elliptic BCm Inozemtsev model. The location of the 
real singularities are completely the same as that in Case 3a. 



7.7 Case 5: A{z) = z(z^ + 2(1 - 2m)z + l)/2 
Change of variable: z — (1 + cn g)/(l — caq). 

Potential: 
V 



+ 





-1) 


4 






m' 


2m 





M 



^ 1 a-{a--l) ^ 1 

f:^ l-cno,: 4 f-tl 



■^-^^ cnQi 4 S 1 + cn Qi 

£- - ^ (A^ - 1 + 2(M - l)c) (A^ + 1 + 2(M - l)c) 

cng. 



M 



+ (A^ + 2(M-l)c)fE 



+ c(c-l)E 



2 ^ dn^ gi 

1^1 — cn cn qj){rn' + m cn g-j cn g^) 
(cn^j — cnqjY 



i=idn 



+ Vb. (7.21) 



Gauge factor: 



M 
i<j 



1 + cn 1 + cn Qj 



1 — cnqi 1 — cn qj 



" ^ a -cng,\ 5 (-^-1+2(^-1)0) 



n 



X 



/ ^ M 



+ cn qi 



exp -TTT arctan — — ^^^^ ) , (7.22) 

V A;A;'(1 - cng^) ' ' ^ ^ 



where A;' = y/m/ — Vl — "t.. The one-body part of the potential is singular 
at integer multiples of 2K{m), so that a natural choice is = 2K(m). This 
model is also regarded as a kind of Jacobi-elliptic BCm Inozemtsev system 
since this case is obtained by a complex projective transformation from Case 
3(m 7^ 0) or Case 4 (cf. Section 9). 



8 Difficulty in Type B Generalization 



Before concluding the article, we will give brief comments on the several- 
variable generalization of the type B monomial space and its difficulty. The 
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type B monomial space of a single variable is defined by 

vi?) = (l,^,...,^-^^). (8.1) 

The general quasi- solvable operators preserving the space (8.1) and the cor- 
responding A/'-fold supersymmetry are investigated in Ref. [22] . It is pointed 
out in Ref. [18] that the type B space (8.1) is realized from the type C space 
(3.6) by the following formal substitution 

jV2 = l, A=M + 1=A^, (8.2) 

although the characteristic features of the type B and C quasi-solvable opera- 
tors are quite different from each other. Hence, a natural generalization of the 
type B space is achieved by the same substitution as Eq. (8.2) in the type C 
space of several variables (3.19), which results in 

= vif2i;M © (<) ■ (8.3) 

Wc easily sec that the above space (8.3) indeed reduces to Eq. (8.1) when M = 
1. In contrast to the single- variable case, however, there is no second-order 
operator preserving the space (8.3) which raises the degree of the monomials by 
two; such an operator must delete simultaneously all the elements cr"^ . . . 
with rij — M — ?> and — 2 as well as the element a-^, but it cannot be 
achieved by any at most second-order operator. In other words, there is no 
several- variable counterpart of an operator which reduces (when M = 1) to 
the degree two second-order operator J_|__|_ preserving the single- variable type 
B space (8.1): 

Needless to say, there are many other possibilities for the generalizations in- 
stead of Eq. (8.3). However, this kind of difficulty may not be the sole obstacle 
and, as far as we have investigated so far, no natural generalization of the 
single-variable type B has been found. 



9 Discussion and Sumniciry 

In this article, we have developed a systematic procedure to construct a fam- 
ily of quasi-solvable operators in several variables leaving two type A spaces 
separately and thus admitting two different polynomial type solutions. It is 
evident that this family is regarded as a subclass of the operators which leave 
a single type A space (3.10) invariant and thus ensure only one polynomial 
type solution. The latter operators are exactly those constructed from the 
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Table 2 

Correspondence between type A and type C canonical forms 



Type 


A 


GL[Z, (L) 


Type C 




Canonical Form 




Case 


T 

1 


1/2 






II 


2z 


Trivial 


1 


III 


2^2 


Trivial 


2 


IV 


2{z^ - 1) 


a = -2/3 = -26 


3(m = 0) 






= -V2e''^/^ 7 = 




V 


2z^ - g2z/2 - 53/2 


P = ei,j = 0,6 = l 


3(m / 0), 4, 5 



In this Table, 52,53 £ C satisfy g2 — 27gl 7^ 0, and ei G C is one of the three different 
roots a {i = 1, 2, 3) of 4/33 - 52/3 - 53 = 0. 

full 5l(M + 1) generators. Therefore, each Hamiltonian of the five cases clas- 
sified here must fall within one of the classes of the sl{M + 1) Lie-algebraic 
Hamiltonians classified in Ref. [13]. To see the correspondence, it is conve- 
nient to extend the free parameters in the Hamiltonians to be complex-valued 
and to consider the GL[2, C) of complex linear fractional transformations in 
Eq. (5.1). 

In Table 2, the first and second columns show the classification scheme and 
the canonical form of the polynomials A{z) under GL{2,C), respectively, for 
the sl{M + 1) Lie-algebraic Hamiltonians preserving one type A space, the 
third column shows the concrete set of the parameters of the GL{2, C) trans- 
formation which renders the each canonical form of the type A to one of the 
canonical form of the type C, and the fourth column shows the corresponding 
type C case classified in this article. For example, we can read from Table 2 
that Case 3 with m = of the type C is equivalent to Case IV of the type A 
under the GL(2,C) transformation while all of Case 3(m 7^ 0), 4, and 5 are 
to Case V. 

Combining the results here with those in Ref. [13], we finally obtain Table 3. 
In Table 3, the number 1 or 2 indicates the number of different families of 
polynomial type solutions available for each model while the symbol x means 
that the corresponding model is only quasi-solvable but not solvable. It is in- 
teresting to note that quantum (quasi-)solvability favors the integrable models 
associated with the BC root system over those associated with the A root sys- 
tem. In the construction of the type C models, on the other hand, we have not 
assumed a priori the invariance under the Weyl group of the BC root system. 
In this respect, we have not appreciated whether this unbalance comes from 
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Tabic 3 

Classification of the (quasi-) solvable quantum many-body systems 



/^l /OKI \ 

• Class (• bubclassj 


Cjuasi-solvable 


C 1 1,1 

solvable 


• Rational A Inozemtsev 




X 


• Rational A Calogero-Sutnerland 




1 


Tl J.* 1 7~>/^ T J. 

• Rational BC inozemtsev 




X 


• Rational BC Calogero-butnerland 




2 


• Hyp.(irig.j A inozemtsev 




X 


•Hyp. (Trig.) A Calogero-butnerland 




1 


+ external Morse potential 






• Hyp. (Trig.) BC Inozemtsev 




X 


•Hyp. (Trig.) BC Calogero-Sutherland 




2 


• Elliptic BC Inozemtsev 




X 


• Elliptic BC Calogero-Sutherland 




X 


• Jacobi-elliptic BC Inozemtsev 


2 


X 



a deeper mathematical structure. 

As was discussed in the previous paper [18], one of the two different famihes 
of polynomial type solutions in most cases does not satisfy the normalizability 
on a space where the Hamiltonian can be naturally defined and thus does not 
provide a physical state. Even in the case where both of them can satisfy the 
normalizability on this space, one of them is still weakly singular at = and 
thus it is natural to discard it as unphysical. As far as the author's knowledge, 
however, there has been no general principle in quantum mechanics which en- 
forces us to do it ^ . Indeed, some recent investigations have rather indicated 
the significance of the careful treatment of this kind of solution, cf. Refs. [24,25] 
and references cited therein. In particular, it was shown in Ref. [24] that by 
employing an ingenious boundary condition at the singularity of the poten- 
tial proper linear combinations of the two linearly independent normalizable 
solutions, which still have a weak singularity, may be well-defined as physical 
eigenstates of the Hamiltonian. Therefore, the possibility of realizing such a 
novel physical state would deserve further investigations both theoretical and 
experimental. 

Although we have mainly concentrated on the construction of quantum many- 
body Hamiltonians in this article, the procedure is quite general and may 



One of the most acceptable postulates could be that in Ref. [23] . 
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have wider applicability. We may construct, for instance, a higher-order linear 
differential operator which admits two or more linearly independent series or 
polynomial type solutions by a slight modification of the present approach. 

Finally, it should be pointed out that the manipulation employed in Section 5 
is justified only after we recognize that the most general type C quasi-solvable 
operator preserves the two subspaces in Eq. (3.19) separately. In other words, 
an analogous condition to Eq. (5.9) in general provides only a sufficient condi- 
tion for an operator to preserve the type C space. In fact, we have found that, 
without the restriction on the value of the parameter A given by Eq. (3.7), 
the most general type C operator no longer preserves the two subspaces sepa- 
rately. The restriction (3.7) was originally introduced [18] in order to prevent 
the type C space (3.6) of a single variable from reducing to the type A space 
(3.9). In the several- variable case, however, it is easily seen that the general- 
ized type C space (3.19) does not reduce to the type A space (3.10) of several 
variables even if the parameter A takes a forbidden value in Eq. (3.7). There- 
fore, there exists the cases, which we shall refer to as irregular cases, where 
the monomial space is of the type C and different from the type A but nev- 
ertheless the general quasi-solvable operator for it does not preserve the two 
subspaces separately. These irregular cases need separate treatment from the 
regular ones presented in this article and will be reported elsewhere. 
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A Formulas 



In this appendix, we summarize the formulas connecting the differential oper- 
ators in terms of Zi with those in terms of 0"^. For the details of the derivation, 
see Ref. [13]. 
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